% nphant
% numerical phantoms

load ../yase/KTR1107_19.mat % calibrated k-trajectory
% characterize imaging experiment
pseq.sw = 103492.884864;    % sampling frequency
pseq.fovcm = 12.8;          % field of view in cm?
pseq.dcoff = 0+i*0;       % DC offset in FID
pseq.offr = 0;             % off-resonance frequency -- units?
pseq.startx = 160;          % ignore this many points at start of FID
pseq.endx = 10;            % ignore this many points at end of FID

% set up image model structure
reso = 64;                  % image model side dimension
wfi = 1.0;                  % relative scaling of exponential imaginary
wfr = 1.0;                  % relative scaling of exponential real
fidfac = 2;
model = model_setup(reso,pseq.fovcm,wfr+1i*wfi);

Ta=.07;               % duration of acquisition in ms
delt=1.0/pseq.sw;           % defined from sw in procpar
NNN=floor(Ta/delt);         % number of acquisition points
NNN=NNN-(pseq.startx+pseq.endx);
tim=(0:NNN-1)'*delt;
N1=3*111;
recon.NLIST=[(N1:N1:13*N1).';NNN];
recon.FLIST=linspace(0.09,0.04,size(recon.NLIST,1)).';
recon.NLENGTHS = size(recon.NLIST,1);  % # of outer loops in PLCG
